System and methods for controlling virtual paths within a network based on entropy rate function

ABSTRACT

A general dimensioning method and system for allocating limited transmission resources to various virtual paths defined on top of a physical network. A two-level hierarchical structure is defined with a layer of one or more virtual paths on top of a layer of physical network elements. Traffic demand is specified for each virtual path and the Entropy Rate Function is used as a blocking measure. The loads on the various links are balanced by equalizing blocking probabilities and the optimal allocation of network physical resources is determined.

CROSS REFERENCES TO RELATED APPLICATIONS

PRIORITY UNDER 35 U.S.C. §119(e) & 37 C.F.R. §1.78(a) (1) This nonprovisional application claims priority based upon a prior U.S. Provisional Patent Application entitled “Efficient Dimensioning Methods For Broadband ATM Networks With General Type of Traffic”, Ser. No. 60/001,169 filed Jul. 14, 1995, in the names of Wlodek Holender and Szabolcs Malomsoky, that is also assigned to the same assignee as the present invention.

This nonprovisional U.S. Patent Application is a continuation of U.S. patent application Ser. No. 08/514,480 entitled “System and Method for Optimal Virtual Path Capacity Dimensioning with Broadband Traffic” filed on Aug. 11, 1995 now U.S. Pat. No. 5,872,918 and contains subject matter related to Nonprovisional U.S. patent application Ser. No. 08/514,235 entitled “System and Method For Optimal Logical Network Capacity Dimensioning With Broadband Traffic”, filed on Aug. 11, 1995 and issued on Jun. 9, 1998 as U.S. Pat. No. 5,764,740 and Nonprovisional U.S. patent application Ser. No. 08/513,723, entitled “System and Method For Adaptive Routing On A Virtual Path Broadband Network”, filed on Aug. 11, 1995, and issued on Mar. 10, 1998 as U.S. Pat. No. 5,727,051. These Nonprovisional U.S. patent applications and the disclosures therein are hereby incorporated by reference herein.

BACKGROUND OF THE INVENTION

1. Technical Field of the Invention

This invention relates to a system and method for the efficient dimensioning of a telecommunications network, and more particularly, to a technique for dimensioning defined virtual paths on a constrained physical network using the Entropy Rate Function as a blocking measure.

2. Description of Related Art

Telephone instruments and other communications devices located in a localized geographic area are conventionally connected with one another by means of switching equipment referred to as local exchanges. Local exchanges, in turn, are interconnected by trunk exchanges. Telephone/data instruments located in geographic areas separated from one another and connected to different local exchanges communicate with one another via a complex grouping of local and trunk exchanges linked together into a telecommunications network. A telecommunication network thus comprises a plurality of interconnected network elements, such as local exchanges, trunk exchanges, mobile radio exchanges, long distance exchanges and combinations thereof. At each network level, traffic from one network element, such as an exchange, to another network element can take various routes that pass through different exchanges.

Efficient network traffic management of the communication facilities within a network requires that a sufficient number of circuits be available to handle the traffic requirements of each destination without exaggerated congestion on last-choice traffic routes. It also requires that the network congestion be as even as possible on all last-choice routes and that unused capacity within the routes which are actually provided be minimized in order to insure efficient utilization of resources. In addition, the telephone company which operates a network has a limited budget and consequently must get as much efficiency as possible out of the existing resources in each network.

In the past, traffic management within a communications network has included procedures for periodically surveying the traffic patterns within the network and changing the configuration of circuits and routes to increase the traffic handling efficiently. In addition, more routes and circuits may be added to a network in anticipation of high call densities to a particular location or a particular region and to local events within that region. Conventional network management systems are also capable of changing the relative distribution of traffic loads between selected routes within the network in order to efficiently maximize the utilization of the network without excessive degradation in the quality of service (QoS). However, conventional traffic network management systems and procedures have generally strived to increase the availability of circuits and routes within a network to handle individual traffic demands rather than to reconfigure a network by redimensioning the routes and circuits within it at a higher level of abstraction so as to maximize the overall efficiency of the network.

In addition to the inherent need to maximize the efficiency of managing physical network resources within a telecommunications system, the growth in information and communications technologies in recent years has created a wealth of new economic opportunities and managerial challenges. Vendors providing telecommunications services are continually faced with new customer demands. Providing ordinary voice communications through a telecommunications network is no longer enough. Today's users want the ability to transmit not just voice signals, but also data, audio, video and multimedia signals in both real time as well as through packet switching networks. Asynchronous Transfer Mode (ATM) technology is acquiring increasing significance due to its enhanced ability to provide broadband telecommunications facilities.

A principal feature of ATM technology is its flexibility in the utilization of network resources. One possible approach to exploit this flexibility is to partition the physical network resources, either completely or partially, into logically defined resources. For example, by partitioning a physical network into a plurality of virtual networks, the operation, maintenance and management of a physical network can be substantially simplified. By this procedure, the task of analyzing and solving the traffic routing problem for each call over an entire physical network can be reduced to the considerably simpler routing problem on each virtual network. Each of the individual virtual networks have less complexity than the entire physical network as a whole, making the solution of the traffic routing problem easier.

Partitioning of physical resources can also be necessitated by the existence of different traffic types, different classes of service or varying traffic demands within leased networks. Instead of providing a separate physical network for each customer application, a service provider can set up a number of virtual networks defined on top of a single common ATM physical infrastructure.

This new flexibility in the configuration of virtual network structures demands efficient dimensioning tools, methods and algorithms. Since the nature of the telecommunications services that will be offered in the future are hard to predict, any dimensioning method that is used to manage the configuration of virtual network resources must be capable of dealing with all types of broadband traffic. If the configuration of a virtual network has to be frequently revised in order to adapt to changes in the pattern of offered traffic, then the computational efficiency of the network redimensioning and configuration control system also needs to be high. An algorithm selected for performing network redimensioning and reconfiguring a network should perform its calculations in a time period that is much shorter than the time duration of each virtual network.

Many existing dimensioning techniques are excluded by the requirement that the dimensioning technique be capable of modeling a general traffic distribution. Most commonly used dimensioning methods are not capable of handling general traffic models because of limitations arising from the use of the Erlang blocking measure. The method and system of the present invention overcomes these disadvantages.

SUMMARY OF THE INVENTION

It is therefore an object of this invention for a user to be able to dimension a network using a general traffic distribution model. It is a further object of this invention to be able to perform the dimensioning calculations using a minimum of computational resources. Another object of the present invention is to implement a dimensioning algorithm that has low computational complexity, and consequently, higher speed. Yet another object of this invention is to be able to perform the redimensioning calculations in as short a period of time as possible. It is a further object of this invention to increase the efficiency of the virtual network dimensioning and the virtual path dimensioning processes.

Given a physical network comprising a plurality of physical links, where each physical link has a prespecified transmission capacity, the system and the method of the present invention illustrates a dimensioning technique that supports a general traffic model. The dimensioning task is treated as a load balancing problem over the various physical links. The optimal solution to the virtual path dimensioning problem corresponds to that choice of allocated capacities over the various virtual paths wherein the blocking on each of the virtual paths is made as uniform as possible over the various links.

In one aspect, the present invention includes a method for efficiently dimensioning a telecommunications network having a plurality of physical links that interconnect a plurality of exchanges or nodes. A plurality of physical links are related to one or more virtual paths. Each of the virtual paths provide an individually switchable connection between a pair of exchanges or nodes in the telecommunications network. Offered traffic is specified for each of the virtual paths and a transmission capacity constraint is set for each physical link of the telecommunications network. The relationship between offered traffic and other computational parameters is modeled on the telecommunications network using an entropy-blocking measure and capacities are allocated to the plurality of virtual paths subject to the transmission capacity constraints for the various physical links such that the blocking probabilities on the various virtual paths are made as uniform as possible within a preselected error bound.

BRIEF DESCRIPTION OF THE DRAWINGS

A more complete understanding of the method and system of the present invention may be obtained by reference to the detailed description of the preferred embodiments that follows, taken in conjunction with the accompanying drawings, wherein:

FIG. 1 is a block diagram of an illustrative telecommunications network within which virtual path dimensioning may be effected;

FIG. 2 is a block diagram illustrating exemplary ATM cell structure;

FIG. 3 is a block diagram illustrating a number of interconnected virtual paths and virtual channels within an ATM network;

FIG. 4 is a block diagram illustrating the cross-connection and switching of virtual paths and virtual channels within an ATM network;

FIG. 5 is a diagram illustrating the CCITT B-ISDN reference model showing the variety of service classes supported and the layers of the standard;

FIG. 6 is a diagram showing an illustrative ATM network providing virtual leased line (VLL) service;

FIG. 7 is a diagram illustrating a multi-layered SDH-based transport network which includes ATM cross-connects;

FIG. 8 is a diagram illustrating the topological relationship between virtual paths and physical links;

FIG. 9 is a flow chart depicting the various steps of an illustrative Push Down algorithm for dimensioning virtual paths defined on a physical network;

FIG. 10 is a graph illustrating the characteristics of the Entropy Rate Function related to the present invention;

FIG. 11 is a graph illustrating the characteristics of the shift parameter that is related to the Entropy Rate Function illustrated in FIG. 10; and

FIG. 12 is a flow chart depicting the various steps of a Virtual Path dimensioning algorithm using the Entropy Rate Function that is related to the present invention.

DESCRIPTION OF THE PREFERRED EMBODIMENT

Congestion Control Within Conventional Networks

Turning first to FIG. 1, there is shown an illustrative schematic diagram of a conventional public telecommunications network including a plurality of local exchanges 21 to 26, each of which have a plurality of local subscribers connected thereto and represented by telephone instruments 27. Two of the local exchanges 21 and 24 are represented as having remote subscriber multiplex stages 28 and 29 associated therewith which, in turn, have local customers 27 connected thereto. The network of FIG. 1 also includes a plurality of trunking exchanges 31 to 34 which serve primarily to interconnect various local exchanges with one another and to provide routes between various parts of the network. Trunk exchange 31 is shown connected to a mobile exchange 35 which includes a pair of illustrative base stations 36 and 37 serving a plurality of mobile radio telephone subscribers represented at 38. In addition, other telecommunications services such as databases and intelligent networks may also be connected to various ones of the exchanges shown. Between each of the exchanges 21 to 35 in the network, there are shown a plurality of communication paths 30, each of which may comprise a plurality of communication circuits, including cables, optical links or radio links for carrying voice and/or data communication between the various exchanges within the network.

The network of FIG. 1 also includes a network control system 40 which is connected to each of the exchanges 21 to 35 within the network by means of communication links 41 (represented by dotted lines) for the transmission of control signals to each exchange and for the reception of traffic data from each exchange. The network control system 40 issues commands to dynamically reconfigure the communication paths within the various traffic routes of the network as well as to control the alarm systems within the exchanges of the network in order to fine tune the alleviation of congestion conditions within the network.

ATM System Concepts

Numerous changes are currently taking place within public telecommunications transport networks which are being implemented. One of the major goals for operators of public telecommunications networks has been to deploy a single type of technology to handle the transport and switching of all types of telecommunications services within a common infrastructure. One such technology is the Asynchronous Transfer Mode (ATM) technology.

ATM is currently being implemented in an attempt to satisfy these needs by creating a bearer telecommunications network which has substantial “bandwidth granularity” and which is capable of coping with very high bandwidth connections. The term “bandwidth granularity” refers to a characteristic of a network in which the bandwidth required by a call can vary freely throughout the duration of the call.

The use of ATM technology in the public telecommunications network provides the capabilities of common switching and transport for related services, increased bandwidth granularity, support of variable-bit-rate services, and support of multimedia services. Because of these features, ATM has been chosen by the International Telegraph and Telephone Consultative Committee (CCITT) as the core technology for broadband ISDN (B-ISDN) services. This is despite the disadvantages of ATM, including transit delays for low speed isochronous services, added complexity within a network, and the introduction of new performance parameters (such as cell-loss and congestion), with which the system of the present invention deals, as will be further set forth below.

An ATM network may be implemented using either plesiochronous digital hierarchy (PDH) or synchronous digital hierarchy (SDH), or both. Moreover, pure ATM may be used as the bearer for a network whenever the limitations arising from multiple conversions between ATM and STM (synchronous transfer mode) and the resultant performance degradations can be dealt with.

The ATM cell structure shown in FIG. 2 is at the heart of ATM technology. An ATM cell has a fixed length of 53 bytes, or octets, divided into a 5-octet header and a 48-octet information field (also known as the payload) . The ATM cell header is structured as a number field and one of its main functions is to assist in routing the ATM cell from the point of origin to the point of destination through one or more switching nodes. The information held in each ATM cell is kept relatively small in order to reduce the size of the internal buffers in the switching nodes and to limit the queuing delays in those buffers. ATM operates in a connection-oriented mode. This is important from a modeling viewpoint since it makes it possible to use the results of well-established circuit-switched mathematical models to optimize the allocation and control of network resources.

The principal function of the ATM cell header is the identification of the virtual connection. Routing information within the ATM cell is contained within two fields: a virtual path identifier (VPI), which determines which virtual path the ATM cell belongs to, and a virtual channel identifier (VCI), which determines which virtual channel in the virtual path the cell belongs to.

A virtual channel is a dynamically allocable end-to-end connection. Optical transmission links are capable of transporting hundreds of megabits per second, whereas virtual channels may fill only a few kilobits per second of a link. Thus, a large number of simultaneous virtual channels can be supported on a single transmission link.

A virtual path, on the other hand, is a semi-permanent connection between endpoints. Each of virtual paths can transport a large number of simultaneously-connected virtual channels. Since a large group of virtual channels are handled and switched together as a single unit, the total processing requirements of a virtual path are less than that of a virtual circuit, and consequently there is faster processing per (virtual) circuit, resulting in a significantly more efficient use of network resources. The network management of virtual paths is relatively simple and efficient.

As illustrated in FIG. 2, the ATM cell header is slightly different at the user-network interface (UNI) compared with the network-node interface (NNI). The UNI contains four bits for generic flow control (GFC) and is used to ensure fair and efficient use of available capacity between a terminal and the network. A payload type indicator (PTI) field is used to indicate whether an ATM cell contains user information or special network information, e.g., for maintenance purposes. A cell loss priority (CLP) field encodes a two-level priority and is used when it becomes necessary to discard cells because of network conditions. The header information is protected by a check sum contained within the header error control (HEC) field.

The use of ATM cells permits the information transfer rate to adapt to the actual service requirements. Depending upon the capacity required, the number of cells per unit of time can be increased up to the transmission bit-rate limit of the physical medium used to carry the data. In addition to data cells there are also cells for signaling and maintenance and idle cells. Signaling cells are used between an end user in the network, or between nodes in the network and their function is to setup a service, e.g., a connection. Maintenance cells provide supervision of the ATM layer while idle cells are used to fill the transmission capacity up to the rate of the transmission medium.

Referring to FIG. 3, there is shown a block diagram illustrating the switching and cross-connection of virtual channels and virtual paths within an ATM link. From the viewpoint of a switch designer, “VP switching” refers to the switching of an ATM cell using only the upper part of the identifier field, that is, the shorter field (VPI) . In contrast, in “VP/VC switching” the entire identified field is used (both VPI and VCI). A VP/VC path consist of a plurality of interconnected VP/VC lengths. Switching and cross-connection can be performed at either the VP or the VC level. The virtual path identifier (VPI) and the virtual channel identifier (VCI) define a two-tier handling and routing structure within the ATM circuitry. From the network architectural standpoint, a virtual path (VP) is a bundle of individual connections, a type of “highway” in the route map of an ATM network. One important task in network management is to allocate the right amount of transmission capacity to each such highway (i.e., a virtual path) in order to optimize network performance. This optimization task is the objective of bandwidth management or virtual path dimensioning techniques and is the subject matter of one aspect of the present invention as further discussed below.

Referring next to FIG. 4, there are illustrated the concepts of virtual path and virtual channel cross-connection and switching. The virtual path identifier (VPI) and virtual channel identifier (VCI) values are only valid for a specific link. In each cross-connect/switch, new VPI/VCI values are assigned to the cell with the combination of physical port and VPI/VCI values providing the identification for the ATM cell. The routing of an exemplary ATM cell is then performed, with the aid of translation tables such as that illustrated in TABLE 1, as follows:

TABLE 1 CONNECTED-FROM CONNECTED-TO PORT VPI VCI PORT VPI VCI A 1 — C 10 — A 2 — D 6 — B 5 3 C 7 4 B 5 2 D 3 5

An ATM cell is the basic multiplexing unit within an ATM transport system, with each cell or information unit containing its own connection and routing information. This feature enables direct multiplexing or demultiplexing of service channels wherein each channel may carry different bit-rates. Each ATM cell is identified and routed by information contained in the header within the virtual path identifier (VPI) and virtual channel identifier (VCI) fields. As mentioned above, a virtual path (VP) is a bundle of multiplexed circuits between two termination points, e.g., switching systems, Local Area Network (LAN) gateways, or private network gateways. A VP provides a direct logical link between virtual path terminations, with the VPI value identifying the particular virtual path.

As also mentioned above, the virtual path concept used within ATM technology allows multiple virtual channels (VCs) to be handled as a single unit. Virtual channels with common properties, e.g., the same quality of service (QoS), can be grouped together in bundles that can be transported, processed and managed as one unit. This flexible bundling simplifies the operation and maintenance of an ATM system.

Both virtual paths and virtual channels can be used to provide semi-permanent paths within the ATM network. Routes are established and released from an operation support system by the setting of “path connect tables” in the cross-connect equipment or in the multiplexers along a path. Virtual channels can also be used for on-demand switching with connections being established by signaling either between a user and the network or within the network.

One important characteristic of ATM technology relates to its protocol architecture and is built around the so-called “core-and-edge” principle. The protocol functions specific to the information type being transported, such as retransmissions, flow control, and delay equalization, are performed in terminals at the “edges” of the ATM network. This leaves an efficient, service-independent “core” network, that only includes simple cell-transport and switching functions. Within the ATM nodes in this core, there are no error checks on the information field nor are there any flow controls. The cell information is simply read, the HEC is then used to correct single-bit errors that might affect the address and the cell is then switched towards its destination.

An ATM adaptation layer (AAL) is used at the edge of the network to enhance the services provided. As shown in FIG. 5, the CCITT reference model for B-ISDN services envisages that the AAL includes service dependent functions. As depicted in FIG. 5, there are three layers in the ATM standard. The first layer is the physical layer defining the physical interfaces and framing protocols. The second ATM layer is independent of the physical medium chosen and defines cell structure, provides multiplexing and demultiplexing and VPI/VCI translation to control the flow of cells within the logical network. The third layer is the AAL which provides the important adaptation between the service and the ATM layer thereby allowing service-independent ATM transport. The AAL performs mapping between the original service format and the information field of an ATM cell. Exemplary functions provided by the AAL include variable-length packet delineation, sequence numbering, clock recovery and performance monitoring.

Deployment of ATM in Telecommunications Networks

One use of ATM technology can be used within customer premises to support high speed data communications in and between customer local area networks. In addition, ATM can be used as an infrastructural resource that is common to all services within a customer premises network, including voice and video communications, data transfers and multimedia applications.

An exemplary service for which ATM nodes are introduced into a public telecommunications network is to provide virtual leased line (VLL) service. VLL service is based upon a virtual path concept and allows line capacity to be directly tailored to customer needs and easily changed without modifying the interface structure. A large number of logical connections can be offered to a user through user-network interfaces (UNIs). In addition, a custom tailored quality of service can also be offered to a customer, matching the needs of the user. Thus, multiple classes of service, quality of service classes and performance parameters can be selected. For example, voice services require low transmission delays but can tolerate high bit-errors, while data communications, on the other hand, are more tolerant of network delays but are sensitive to bit-errors. Thus, the quality of service level of a particular application can be contractually agreed to between a service provider and a customer and audited manually or automatically to ensure compliance.

As shown in FIG. 6, there is an exemplary virtual channel based VLL service implemented within a ATM network. Network terminals A to E are each coupled through flow enforcement nodes 601 to 605, respectively, to ATM cross-connect nodes 611 to 614. The ATM network consist of a plurality of ATM cross-connects 611 to 614 which can provide routing both at the virtual path as well as at the virtual channel level. The flow enforcement functions 601 to 605 are located at the edge of the ATM network to protect the network against potential overloads. This function ensures that no connection violates the conditions agreed-to when the connections are setup. Additional services can be implemented by adding services to one or more of the cross-connect nodes 611 to 614. Within the network of FIG. 6, an exemplary virtual path is illustrated by the wavy line 621 between terminal C and D. A first virtual connection between terminals A and B is illustrated by the dashed line 631 while a second virtual connection is illustrated by the dotted line between terminals C and E, as illustrated by the dotted line 632.

In addition to the virtual leased line network shown in FIG. 6, other services, such as SMDS/CBDS and frame relay, can easily be added depending upon demand by connecting servers to the ATM nodes within the network. In residential areas, ATM technology can be used to provide new and enhanced entertainment services such as on-demand video to the end user. The flexibility of an ATM network makes it possible to support a multitude of services, such as long distance education, home shopping, and games.

FIG. 7 illustrates an ATM network which has been overlaid upon a SDH-based layered transport network. The layers include a customer premises network layer 701, a local transport network layer 702, a regional transport network layer 703 and a national transport network layer 704. A plurality of ATM business network nodes 711 to 714, control the flow of data from the customer premises terminals 715 and LANs 716 into respective ones of a plurality of add-drop multiplexers (ADM) 721 serving SDH cross-connect nodes 722 within the local transport network 705. The local cross-connect nodes 722 are in turn coupled through regional cross-connect nodes 731 in the regional transport network, two of which are coupled by add-drop multiplexers 732. Within the local transport network layer 702, a pair of ATM access nodes 723, and SDH rings, comprising the add-drop multiplexers 721, serve the cross-connects 722 and are used for subscriber access with a capacity of up to a full STM-1 (155 megabits per second), the standardized access rate for B-ISDN services.

Existing traffic such as Plain Old Telephone Service (POTS) can also be carried on this ring network, with remote multiplexers and other access nodes providing the final local-loop connection. The ATM access nodes 723 are shared for access to different services from one location and can include both voice and data using different VP/VCs. In the ATM access nodes 723, ATM traffic is concentrated to make more efficient use of the transport capacity.

The size of a ATM access node can vary, depending upon the capacity required, from a small multiplexer to a large cross-connect. In the regional transport layer 703, ATM cross-connects 733 are used to route traffic between local areas. In the national transport network layer 704, illustrated in FIG. 7, ATM is not visible. In an ATM overlay network, in place as illustrated in FIG. 7, services such as frame relay and SMDS/CBDS are easily added. Functionality for B-ISDN can also be added to both access and regional nodes by adding appropriate software and hardware. As also illustrated in FIG. 7, a network management system 750, such as one operating in accordance with the TMN standard of the CCITT can be implemented to provide the necessary network management functionality to both the SDH and ATM elements of the network.

The management of the ATM network by subsystem 750 may be implemented in accordance with the telecommunications management and operations support (TMOS) family of network management systems provided by Telefonaktiebolaget LM Ericsson, the assignee of the present application. Such network management may include various functionalities such as routing algorithms and congestion control implemented in accordance with the teachings of the present invention as set forth in detail below.

Virtual Path Capacity Dimensioning

A useful model in dimensioning a telecommunications network is to treat the dimensioning problem as involving a two-layered structure consisting of a first physical network layer, with discrete topology and specified link capacities and a second virtual path layer with virtual paths and their particular routings. A traffic demand is offered to the virtual paths in this model. In dealing only with the task of dimensioning network capacity, the virtual paths are, ipso facto, already routed. Each virtual path may travel through numerous physical links but will emulate a highway comprising only a single path. Each virtual path will have one characteristic blocking value and one characteristic allocated capacity value with only as many variables in the model as there are virtual paths.

The term “offered traffic” is used to refer to the time varying demand for transmission capacity along each virtual path. The term “traffic demand” is used to indicate the time-averaged value of the offered traffic for each link. When the characteristics of the traffic on an ATM network can be modeled by a single-parameter Poissonian distribution, the traffic is termed as homogenous single-class traffic. When the offered traffic is non-homogeneous, it is usually modeled using a multi-class Poissonian distribution.

The offered traffic can also be modeled by a normal distribution. This is referred to as normal traffic. Finally, network dimensioning can also be based upon actual traffic as determined by measurements.

The transmission requests of a multitude of users can be combined into an aggregated traffic stream. For example, several users could be sending messages from Dallas to Stockholm contemporaneously. Instead of managing these multiple transmissions individually, it is more efficient to combine and transmit them as a group over a broadband trunk line. As discussed above, a virtual channel is a dynamically allocable end-to-end connection. A virtual path is a logical construct that allows multiple virtual channels to be handled and switched together as a single unit. This unified switching reduces the processing requirements overall and also speeds up the transmissions. Since the management of a virtual path is both simpler and more efficient than the management of virtual channels or individual physical circuits, a significantly improved utilization of network resources can be obtained by this technique.

The Virtual Path Dimensioning Model

The basic model that is considered is that of a connection-oriented network that is operating under a fixed routing. If we define the physical network as comprising a set J of arbitrarily connected links, then each virtual path (VP) or route r is an ordered list whose elements are a subset of J. The relationship between the virtual paths and the physical links can be defined in terms of a routing matrix, X, whose elements are: $\begin{matrix} {\chi_{j,r} = \left\{ \begin{matrix} 1 & {{{if}\quad j} \in r} \\ 0 & {otherwise} \end{matrix} \right.} & \left( {{EQN}.\quad 1} \right) \end{matrix}$

FIG. 8 illustrates the topological relationship between virtual paths and physical links. In FIG. 8, virtual path, VP₁, consists of physical links P₁, and P₂, and virtual path, VP₂, consists of physical links P₂ and P₃.

The quantitative relationship between the capacities allocated to various VPs and the corresponding capacities allocated on physical links is given in matrix form as:

χC _(VP) =C _(phys)  (EQN. 2)

where χ is a routing matrix of the type defined in Equation (EQN. 1) above, C_(VP) is the virtual path capacity vector and C_(phys) is the physical link capacity vector.

The physical link capacity vector, C_(phys), representing the capacities that have been allocated to the various VPs, cannot exceed the available physical capacities on any of the physical links. This limitation can be expressed by the simple constraint relation:

χC _(VP) =C _(phys) ≦{overscore (C)} _(phys)  (EQN.3)

where {overscore (C)}_(phys) is a vector describing the available physical capacities on each of the K physical links. Since Equation (EQN. 3) is a vector inequality, it is important to note that corresponding components of both vectors have to satisfy this inequality. For the simple example shown in FIG. 8, the routing matrix χ is: $\begin{matrix} \begin{bmatrix} 1 & 0 \\ 1 & 1 \\ 0 & 1 \end{bmatrix} & \text{(EQN.~~4)} \end{matrix}$

where the relation between virtual paths capacities and corresponding physical link capacities allocated to VPs is: $\begin{matrix} {{\begin{bmatrix} 1 & 0 \\ 1 & 1 \\ 0 & 1 \end{bmatrix} \cdot \begin{bmatrix} C_{1} \\ C_{2} \end{bmatrix}_{VP}} = \begin{bmatrix} C_{1} \\ {C_{1} + C_{2}} \\ C_{2} \end{bmatrix}_{phys}} & \text{(EQN.~~5)} \end{matrix}$

where ${\begin{matrix} C_{1} \\ C_{2} \end{matrix}}_{VP}$

is the virtual path capacity vector and ${\begin{matrix} C_{1} \\ {C_{1} + C_{2}} \\ C_{2} \end{matrix}}_{phys}$

is the allocated vector for the physical links.

The call request process to a given route r can be any stationary process for which we know the fictitious occupancy distribution if the process were to be provided with a resource having infinite capacity that accepted all calls and avoided all blockage. X_(r) denotes the occupancy level of this fictitious infinite capacity resource, and is commonly referred to in the art as “offered traffic”.

The virtual path dimensioning problem is defined in the system and method of the present invention as having two objectives. First, the transmission capacity allocated to each virtual path needs to be optimized so as to minimize a transmission cost function. Second, for every physical link, the capacities allocated to the various virtual paths traversing this link must not exceed the physical transmission capacity constraint of that physical link.

The physical capacities that can be allocated to the various virtual paths can be approximated by a continuous function that assumes any real value in the range [0, {overscore (C)}_(phys)]. Consequently the resource optimization task does not require discrete optimization and all of its attendant complexities. In the dimensioning problem addressed in this application, load sharing between different VPs has not been considered. The offered traffic is assumed to be defined for each virtual path. Further, because the network has fixed routing, the routing of the offered traffic is fixed by the choice of the virtual paths.

The “Push Down” Dimensioning Technique

The virtual path dimensioning task is viewed in the present invention as a load balancing problem in which the “load” is the value of an appropriate chosen blocking measure and where the optimal solution corresponds to that choice of allocated VP capacities wherein the blocking on each of the virtual paths is distributed in as uniform a manner as possible. One way to make the blocking distribution uniform is to measure the divergence in the values of the blocking on the various virtual paths and then to minimize this divergence. This approach can be implemented using any standard minimization algorithm, for example, the well-known simulated annealing technique.

A related approach would be to first identify the virtual path having the highest blocking value and then to minimize the blocking for this virtual path by reallocating capacity from other VPs until the virtual path is no longer the VP with the highest blocking. This formulation corresponds to a min-max optimization problem and can be analytically formulated as described below.

If we denote the blocking on the i^(th) virtual path as B(VP_(i)), then the VP that has the largest blocking is max(B(VP_(i))) where the maximum is taken over all the VPs. The maximum of a blocking measure over the set of virtual paths defines the objective function (also known as the cost function) for the VP dimensioning problem. The goal of the optimization procedure therefore, is to find the minimum of the objective function, which corresponds to:

min max (B(VP _(i)))  (EQN. 6)

where the minimum is defined over all feasible configurations.

Since this technique involves pushing down the highest blocking value among all the VPs considered, an algorithm using this technique to solve an optimization problem is called an “Push Down” algorithm. This algorithm follows from the fact that the uniform blocking distribution corresponds to the best solution of the unconstrained VP dimensioning problem. Accordingly, the best solution is to allocate capacities to each VP such that the blocking on each of the VPs is made equal within an error bound. However, such a solution is not always realizable because of the capacity constraints of the various physical links. The limited capacity of a physical link needs to be shared amongst all VPs traversing that physical link.

FIG. 9 shows the various steps in one illustrative embodiment of a Push Down Algorithm for dimensioning virtual paths defined on a physical network. The dimensioning process starts at 902 with the definition of the connection topology of the various VPs. The various VPs are also assembled into a VP Dimensioning Set. The VPs are then grouped at 903 in order of the physical links that they each traverse. Initial allocations of transmission capacity are then made to each VP at 904. At 905, a target FIG. for reduction in blocking is chosen. To set a target, one needs to first select a blocking measure. In one preferred embodiment of the present invention, the Entropy Rate Function, as detailed in the sections following, is used as a blocking measure. The target value is used to set termination conditions for the dimensioning algorithm.

The blocking on each of the VPs is determined at 906 over each of the physical links. If the various VPs traversing a single physical link do not face the same or similar levels of blocking, then the capacities currently allocated to each of the VPs are revised at 907 so as to equalize the blocking values for the VPs within an error bound. Capacity can be added to VPs by the allocation of unallocated physical capacity or by reallocation of already allocated capacity from a less productive VP to a more productive VP. This capacity readjustment is performed without violating the capacity constraints of any of the physical links.

This process results in the identification at 908 of one or more physical links as being the bottle-necks in this optimization procedure. A physical link on which the VP blockage is the highest and whose blockage is not reducible by capacity reallocation is referred to as a critical link. Each critical link determines the lowest blocking that can be achieved on VPs traversing that physical link. One of the principal tasks of the Push Down algorithm is to identify, at each stage of the optimization procedure, the set of critical links for a given set of virtual paths.

Once a critical link is identified at 908, physical capacities can be reallocated between the various virtual paths traversing this critical link in such a way as to equalize the blocking values for each of the virtual paths. It should be noted that when a physical link is found to be a critical link, ipso facto, it has no unallocated capacity. Consequently, only reallocation of capacity between VPs passing though a critical link is possible after the algorithm reaches this stage of the dimensioning procedure.

Next, at 909, VPs with allocated capacities are removed from the set of all the VPs that still need to be dimensioned. Correspondingly, at 910 the available physical link capacities are reduced by the allocated capacities of the VPs that were removed in the previous step.

The dimensioning problem thus reduces to the optimization problem of minimizing the highest blocking probability for the remaining set of VPs. This permits the use of a recursive re-entrant algorithm to implement this procedure.

The blocking values from the previous step are now used as the initial values in the remaining dimensioning problem. This optimization procedure is recursively repeated at 911 until all the capacities of each of the physical links have been allocated. In summary, this greedy-type algorithm starts with dimensioning the complete set of all VPs and terminates at 912 when the set of virtual paths remaining to be dimensioned becomes a null set.

It should be emphasized that the execution of every dimensioning algorithm of the type detailed herein need not identically follow the order of steps as depicted in FIG. 9. The execution order of some of the steps of the dimensioning algorithm may be different from that shown in FIG. 9 based on implementation details and convergence considerations.

The problem of analytically identifying the critical link in a given set of VPs, has proven to be a difficult task. There are no known techniques for determining the critical link directly from the offered traffic and the physical link capacity constraints. Hence the push down algorithm employs an iterative approach to identify critical links. The algorithm is initialized for all VPs by using a uniform large blocking value for all the VPs. The initial blocking value that is selected has to be large enough so that the sum of the initially allocated values of the VP capacities do not exceed the available physical capacities of the various physical links.

By the slow and uniform reduction in the degree of blocking on the set of all the virtual paths remaining in the optimization procedure, the critical link is identified at each level as being that link which first violates the physical capacity constraints of a traversed physical link.

Dimensioning Using an Entropy Blocking Measure

The speed and efficiency of the above procedure for identifying the critical link at each stage of the dimensioning process is critically dependent upon the complexity of the blocking measure used in the modeling. Traditionally, the Erlang blocking measure (also known as the time congestion blocking formula) has been used to determine the optimal allocation of VP capacities in a network.

The present technique incorporating the use of the Entropy Rate Function as a blocking measure yields superior results to those obtainable by use of the Erlang blocking measure. The use of the Entropy Rate Function permits the modeling of arbitrary traffic distributions, and in most cases this computation can be done much faster compared to computations based upon other blocking measures. It has also been found that the iterative search for the critical link can be substantially improved, a result principally following from the fact that the Entropy Rate Function is a convex function. Prior to a description of a dimensioning algorithm using the Entropy Rate Function, it would be useful to explore the characteristics of the Entropy Rate Function.

The Entropy Rate Function as a Blocking Measure

As noted earlier, the choice of the blocking measure is critical to the Push Down algorithm. A general expression for the blocking measure based on the Entropy Rate Function will be derived next and applied to an exemplary situation where the offered traffic is alternately modeled by single-class and by multi-class Poissonian distributions.

The Entropy Rate Function is known in the art, and has been used to model congestion at the physical link level, see, e.g., J. Y. Hui, A Congestion Measure for Call Admission and Bandwidth Assignment for Multi-Layer Traffic, International Journal of Digital & Analog Cabled Systems (1990), but has not hitherto been used as a blocking measure in solving either the dimensioning or the planning problem at either the virtual path level or at the network level. Additionally, the Entropy Rate Function has been used to define the concept of “effective capacity” of a physical link. It is important to note that the dimensioning technique using the Entropy Rate Function detailed herein is not limited to offered traffic that follows a Poisson distribution and that the system and method works equally well with any type of offered traffic distribution, including that determined by measurements.

Saturation blocking probability can be defined as the probability that the traffic demand exceeds a specified value of transmission capacity. The saturation probability is also called the “tail probability”, because it denotes the probability mass of the tail of the offered traffic distribution. A well known approximation to this tail probability, namely Chernoff's Bound, is derived below.

Let X be an arbitrary distributed random variable, and C be a given value. Then for all positive values of s, the following bound can be shown to exist: $\begin{matrix} {{P\left( {X > C} \right)} = {{{P\left( {^{sX} \geq ^{sC}} \right)} \leq \frac{E\left( ^{sX} \right)}{^{sC}}} = ^{- {({{sC} - {\ln {({E{(^{sX})}})}}})}}}} & \text{(EQN.~~7)} \end{matrix}$

where P(X>C) is the probability that the random variable X takes values greater than C.

The derivation of this bound is based on the Markov inequality. The term, ln(E(e^(sx))), denotes the log moment generating function, and is also called the cumulant function, μ(s). The tightest bound (also known as Chernoff's Bound) is obtained by maximizing the exponent sC−μ(s) with respect to s. This maximum is reached at s=s* where s* is the unique positive solution for the equation, C=μ′(s).

It is easy to show that μ′(s) increases with s, thereby guaranteeing the uniqueness of its root. μ′(s) increases with s because μ″(s)≧0. This in turn is due to the fact that the second derivative is equal to the variance of the (shifted) distribution. This maximizing exponent is denoted by I_(x)(C) and is called the Entropy Rate Function. The Entropy Rate Function is the convex conjugate transform of the cumulant function and can be stated by the following expression: $\begin{matrix} {{{{{- \ln}\quad {P\left( {X \geq C} \right)}}\quad \geq \quad {\sup\limits_{s}\quad \left( {{sC} - {\ln {\quad \quad}{E\left( ^{sX} \right)}}} \right)}} = {I_{x}(C)}}\quad} & \text{(EQN.~~8)} \end{matrix}$

Since lnP(X≦C)=lnP(−X≧−C), and lnP(−X≧−C)=lnP(e^(−sX)≧e^(−sC)), a similar expression for the left tail of the distribution can be derived by the following inequality: $\begin{matrix} {{{{{- \ln}\quad {P\left( {X \leq C} \right)}}\quad \geq \quad {\sup\limits_{s}\quad \left( {{- {sC}} - {\ln {\quad \quad}{E\left( ^{- {sX}} \right)}}} \right)}} = {I_{- x}\left( {- C} \right)}}\quad} & \text{(EQN.~~9)} \end{matrix}$

The relation between the right and the left tails of the Entropy Rate Function distribution can be expressed using the parameter, s as:

I _(−X)(−C(s))=I _(X)(C(−s))  (EQN. 10)

Thus, by changing the sign of the parameter, s, one can switch from the right tail to the left tail of the Entropy Rate Function distribution and vice versa.

The Entropy Rate Function for Homogenous Poissonian Traffic

The use of the Entropy Rate Function to dimension virtual paths when the offered traffic is homogenous is considered first. Homogenous Poisson traffic can be characterized by a bandwidth demand parameter, p, the mean call arrival rate, r, and the mean duration of each call, h, (also referred to as the mean holding time). The traffic demand, ρ, is therefore the product of the mean call arrival rate and the mean holding time, i.e., r*h. The cumulant function for homogenous traffic is described by the relationship:

μ(s)=ρ(e ^(sp)−1)  (EQN. 11)

Consequently, the allocated capacity, C, and the Entropy Rate Function, I, for homogenous traffic are given by:

C=μ′(s)=ρpe ^(sp)  (EQN. 12)

and

I(C)=s*(C)·C−μ(s*(C))  (EQN. 13)

i.e.

I(C(s))=sρpe ^(sp)−ρ(e ^(sp)−1)  (EQN. 14)

By substituting the relation $s = {\frac{1}{p}\left( {{\ln \quad C} - {\ln \quad \left( {\rho \quad p} \right)}} \right)}$

as the solution to Equation (EQN. 12), the Entropy Rate Function for the homogeneous traffic of the type described by Equation (EQN. 14) can also be expressed solely as a function of the allocated capacity, C, as: $\begin{matrix} {{I(C)} = {{\frac{C}{p}\left( {{\ln \quad C} - {\ln \quad \left( {\rho \quad p} \right)}} \right)} - {\rho \quad \left( {\frac{C}{\rho \quad p} - 1} \right)}}} & \text{(EQN.~~15)} \end{matrix}$

FIG. 10 graphs the characteristics of the Entropy Rate Function for different values of offered traffic in the unit bandwidth demand, p. FIG. 11 shows the shift parameter, s, plotted as a function of the allocated capacity, C.

As shown in FIGS. 10 and 11, the Entropy Rate Function has three significant properties. First, it is a convex function that reaches its minimum value of zero at the mean of the distribution, i.e., when C=ρ. Second, the shift parameter, s, goes from negative to positive for increasing values of C at the mean of the distribution, i.e., at C=ρ. As can be seen from FIG. 11, the shift parameter, s, is negative when C<ρ and is positive when C>ρ. Third, the shift parameter, s, is monotonic and increases as a function of the capacity allocated to a virtual path.

The transform parameter, s, can thus be interpreted as a probability distribution shift parameter. When the shift parameter takes negative values, the probability distribution is shifted to the left compared to the probability distribution corresponding to the zero value of the shift parameter. When the shift parameter takes positive values, then the probability value is shifted to the right.

The Entropy Rate Function for Multi-Class Poissonian Traffic

The traffic model can be extended to offered traffic characterized by a multi-class Poisson distribution and the Entropy Rate Function corresponding to such an offered traffic model is derived next.

Replacing the entropy measure for the single-class distribution with the entropy measure for the multi-class distribution induces the difficulty that the Entropy Rate Function can no longer be explicitly expressed in terms of the allocated capacity, C. To avoid this problem, the Entropy Rate Function is expressed in terms of the shift parameter s which is also utilized as a control parameter. By increasing the absolute value of this parameter, allocated capacities can be varied implicitly, thus permitting the entropy measure to be positively incremented.

Let us denote the random traffic of class i (where i takes values from 1 to k) by X_(i). The peak bandwidth demand p_(i), the mean call arrival rate r_(i) and the holding time h_(i) of each class are defined such that ρ_(i)=r_(i)h_(i) and the expectation value of the cumulative value of the random variable X_(i) is: $\begin{matrix} {{E\left( {\sum\limits_{i = 1}^{k}X_{i}} \right)} = {\sum\limits_{i = 1}^{k}{\rho_{i}p_{i}}}} & \left( {{EQN}.\quad 16} \right) \end{matrix}$

The Entropy Rate Function $I_{\sum\limits_{i = 1}^{k}\quad X_{i}}(C)$

for multi-class traffic $\sum\limits_{i = 1}^{k}X_{i}$

is an estimate of the negative logarithm of the saturation probability: $\begin{matrix} {{- {\ln \left( {{\sum\limits_{i = 1}^{k}X_{i}} \geq C} \right)}} > {I_{\sum\limits_{i = 1}^{k}\quad X_{i}}(C)}} & \left( {{EQN}.\quad 17} \right) \end{matrix}$

where k is the number of traffic classes and C is the capacity allocated to the VP carrying this aggregated multi-class traffic.

The following relation expresses a general property of the Entropy Rate Function that is valid for all distributions: $\begin{matrix} {{I_{\sum\limits_{i = 1}^{k}\quad X_{i}}\left( {E\left( {\sum\limits_{i = 1}^{k}X_{i}} \right)} \right)} = {{I_{\sum\limits_{i = 1}^{k}\quad X_{i}}\left( {\sum\limits_{i = 1}^{k}{\rho_{i}p_{i}}} \right)} = 0}} & \left( {{EQN}.\quad 18} \right) \end{matrix}$

It can likewise be shown that $I_{\sum\limits_{i = 1}^{k}\quad X_{i}}(C)$

is convex in C.

The log moment generating function for multi-class Poisson traffic is given by the relation: $\begin{matrix} {{\mu (s)} = {\sum\limits_{i = 1}^{k}{\rho_{i}\left( {^{{sp}_{i}} - 1} \right)}}} & \left( {{EQN}.\quad 19} \right) \end{matrix}$

Since the Entropy Rate Function has the general form: $\begin{matrix} {{I_{X}\left( {C(s)} \right)} = {{\sup\limits_{s}\quad \left( {{sC} - {\mu_{X}(s)}} \right)} = {{s\frac{}{s}{\mu_{X}(s)}} - {\mu_{X}(s)}}}} & \text{(EQN.~~20)} \end{matrix}$

the Entropy Rate Function for multi-class traffic can be expressed in terms of the shift parameter, s, as: $\begin{matrix} {{{I_{\sum\limits_{i = 1}^{k}X_{i}}\left( {C(s)} \right)} = {\sum\limits_{i = 1}^{k}\left( {{s\quad \rho_{i}\quad p_{i}^{{sp}_{i}}} - {\rho_{i}\left( {^{{sp}_{i}} - 1} \right)}} \right)}}\quad} & \text{(EQN.~~21)} \end{matrix}$

where the allocated capacity, C, can further be functionally related to the shift parameter, s, as: $\begin{matrix} {{C(s)} = {{\frac{}{s}{\mu (s)}} = {\sum\limits_{i = 1}^{k}{\rho_{i}\quad p_{i}^{{sp}_{i}}}}}} & \text{(EQN.~~22)} \end{matrix}$

The replacement of the single-class entropy measure with the multi-class entropy measure makes the problem analytically more complicated because the entropy measure can no longer be expressed explicitly in terms of the allocated capacity, C. This complication is due to the fact that the shift parameter, s, cannot be eliminated from Equation (EQN. 22).

However, since Equation (EQN. 21) expresses the Entropy Rate Function in terms of the shift parameter, s, we can vary s instead of varying C. Thus the capacity value can be computed at each iterative step of the algorithm by using Equation (EQN. 22). It should be noted that the zero value of s corresponds to the zero value of the entropy measure. The dimensioning algorithm is initialized by setting s to zero for all VPs.

The Entropy Rate Function for Normally Distributed Traffic

The Entropy Rate Function can also be used with other traffic models than those based on Poissonian offered traffic distributions. Two other important traffic models are discussed below. The first traffic model is based on the Entropy Rate Function for a normal distribution of offered traffic. Corresponding expressions for the Entropy Rate Function are derived for this traffic model. The second traffic model is based on the Entropy Rate Function derived from actual traffic flow measurements rather than on explicit assumptions about the distribution of the offered traffic.

The Entropy Rate Function for normally distributed traffic has been shown, see R. S. Ellis, Entropy, Large Deviations and Statistical Mechanics 39 (Springer-Verlag, 1985), to be defined by the following expression: $\begin{matrix} {{I_{X}(C)} = \frac{\left( {C - m} \right)^{2}}{2\quad \sigma^{2}}} & \left( {{EQN}.\quad 23} \right) \end{matrix}$

where m is the mean and σ is the variance of the normal distribution N(m, σ). Further, $\begin{matrix} {{s(C)} = \frac{C - m}{\sigma^{2}}} & \left( {{EQN}.\quad 24} \right) \end{matrix}$

These two relations yield the following simple expression for the Entropy Rate Function directly in terms of the control parameter, $\begin{matrix} {{I_{X}\left( {C(s)} \right)} = {\frac{\sigma^{2}}{2}\quad s^{2}}} & \left( {{EQN}.\quad 25} \right) \end{matrix}$

Thus in the case of normally distributed traffic, the Entropy Rate Function is shown to be a simple (and convex) quadratic function.

The Entropy Rate Function for Measured Traffic

Since the diversity of services that will be provided in the networks of the future will be considerably greater than those available today, it is very likely that future broadband networks will experience huge variations in the types of traffic offered to the network. Consequently, all models that are based upon assumptions about specific idealized representations of the traffic distribution are likely to be inadequate because of their inherent inflexibility. Using an Entropy Rate Function derived from traffic measurements may provide one solution to this difficult traffic model estimation problem.

The traffic models that were discussed up to now had been defined on the call-level time scale. In contrast, traffic measurement statistics are defined on the time scale of a standard ATM cell. The call-level time scale can be considered to be an approximation of the cell-level time scale. Randomly changing traffic demand during a call can thus be abstracted by one or more parameters describing constant bandwidth demand on the call-level time scale.

It has recently been suggested that the Entropy Rate Function can be estimated from traffic measurements on the cell level. See, e. g., N. G. Duffield et al., Entropy of ATM Traffic Streams: A Tool for Estimating QoS Parameters (Dublin Institute for Advanced Studies, 1994).

Effect of Traffic Overload on the Entropy Rate Function

The interpretation of the Entropy Rate Function as a blocking measure, works well if the average offered traffic on each physical link is lower than the corresponding available physical capacity on that link, i.e., if E(X_(k))≦{overscore (C)}^(k) _(phys). However, this condition may be violated in some realistic overload situations. Consider the following example based on homogeneous Poisson traffic and time-congestion blocking measure (i.e., the Erlang blocking formula).

TABLE 2 lists three values of allocated capacities and corresponding traffic demands that were calculated for a fixed blocking value of 0.03. Note that in the last case, the offered traffic is larger than the allocated capacity, even though the blocking is relatively low.

TABLE 2 ALLOCATED CAPACITY AVERAGE OFFERED TRAFFIC 15.45 10 109.29 100 699.00 700

This example shows that the entropy blocking measure needs to be extended to cover the overload situation when the condition E(X_(k))≦{overscore (C)}_(phys) is violated. Mathematically, such an extension can be accomplished easily. As shown previously, the Entropy Rate Function is a convex function having a minimum value of zero at E(X_(k)) . The left branch of the Entropy Rate Function defines the overload region (see FIGS. 10 and 11). In this region, an increase in the Entropy Rate Function corresponds to a decrease in the allocated capacity and to negative values of the control parameters. By changing the sign of the control parameter, Push Down algorithm based on the Entropy Rate Function can be easily extended to cover the overload region. Such an extension requires only a minor modification of the original dimensioning algorithm.

The problem that remains is largely conceptual in nature, namely how to interpret this extension. The left branch region of the Entropy Rate Function corresponds to the approximation of the left tail of the probability mass:

−lnP(X≦C)≧I _(X)(C(−s))  (EQN. 26)

and can be interpreted as an entropy measure of the utilization of the allocated capacity, C.

Since in the initial state, the resources that were allocated exceeded the capacities of the physical resources, i.e., E(X_(k))≧{overscore (C)}_(phys), consequently, the utilization has to be decreased as it corresponds to an increase in the Entropy Rate Function.

One interpretation of the objectives of the optimization in the overload region can be as follows. To improve the uniformity of the distribution of entropy utilization measures in this region, the largest consumer of resources (i.e., the VP having the lowest entropy) is identified and the utilization of this extreme representative is decreased. Decreasing utilization of the largest consumer corresponds to an increased Entropy Rate Function in the overload region. Thus, this approach corresponds to a maxi min formulation of the optimization problem. It should be noted that the constraint conditions are now applied by starting with values that initially fall in the overload region.

Here again, we can use the correspondence of the uniform distribution of the entropy utilization measures to the best utilization of resources can be used (even if the best utilization is infeasible). Following the same reasoning as for the right branch region, to satisfy the constraint conditions, the entropy utilization measure is made uniform on each critical link. Further, the capacities of all VPs traversing the critical link are allocated in such a way that the entropy utilization values of these VPs are made equal.

As discussed earlier, the objective of the optimization in the right side region of the entropy curve was to increase the capacity allocated to the VP having the highest blocking (i.e., the VP with the least entropy blocking measure). This corresponds to a maxi min formulation of the optimization problem. It should be noted that the optimization objective for the left side region can be transformed to an optimization problem for the right side region by replacing the term “utilization” by the term “blocking” and the term “largest consumer of resources” by the term “VP with the largest blocking” in the earlier formulated optimization objective for the left side region.

These two different formulations of the optimization objective for the left and the right branch of the Entropy Rate Function result in the same optimization procedure. In both cases the Entropy Rate Function needs to be increased. This can be accomplished by increasing the absolute value of the control parameter, s. If the load does not exceed the resources available, then the shift parameter will be positive and the capacities allocated to the various virtual paths can be successively raised until all the available physical resources are allocated. If, on the other hand, the load exceeds the resources available, then the shift parameter is negative. In such cases, the allocated capacities need to be gradually decreased until they come within the physical resource constraints.

VP Dimensioning Algorithm Using the Entropy Rate Function

We can apply the above-detailed characteristics of the Entropy Rate Function to solve the VP dimensioning problem efficiently. As explained earlier, the VP dimensioning problem aims at allocating limited physical network resources amongst a plurality of predefined VPs given an offered traffic distribution. One embodiment of a VP dimensioning algorithm that uses the Entropy Rate Function as a blocking measure is shown in FIG. 12.

The process starts with a series of initialization steps 1202 to 1206. All of the VPs that are to be dimensioned are assembled into a VP dimensioning set at 1202. The transmission capacity constraints for each physical link in the network are specified at 1203. A set of upper limits, for the Entropy Rate Function, I_(MAX), are optionally specified, one per virtual path, at 1204.

It should be noted that steps 1202 to 1204 may be performed in any order depending on implementation considerations. Further, I_(MAX) is specified at 1204 only if there is a realistic likelihood of the offered traffic distribution having a truncated right tail i.e., only if P(X>C) is zero for values of X that are greater than some finite value, C_(MAX). If an offered traffic distribution has a truncated right tail, then it is theoretically possible to size network resources to achieve zero blocking. However, such situations are rare in practice.

Other initialization steps include the selection of large and equal blocking values for each of the Virtual Paths at 1205. As explained elsewhere, there is an inverse relationship between the values of the Entropy Rate Function and the corresponding blocking on a VP. Consequently, large blocking values correspond to small values of the Entropy Rate Function. Using the relations developed earlier, initial capacity allocations for the various VPs are also computed at 1205.

At 1206, these initial capacity allocations are cumulated over each physical link and compared to the prespecified transmission capacity of that physical link. If the initial allocations are such that there is overload on one or more of the physical links, then the computational steps that follow are modified to use formulae based on negative values of the shift parameter, s. This result follows from the fact that the algorithm for the overload situation corresponds to a reflection of the dimensioning problem to negative values of the shift parameter, s.

Another step implicit in the initialization process is to choose an offered traffic model, and if this model is not based on measurements, to derive the corresponding relations for the Entropy Rate Function, I, the allocated capacity, C, and the shift parameter, s. This step is not shown in FIG. 12.

After the initialization steps 1202 to 1206, the dimensioning technique performs recursive steps 1207 to 1215. The recursive technique outlined in FIG. 12 shows a two-level recursion wherein the VP dimensioning algorithm first iteratively allocates capacity to the VPs in the VP dimensioning set, as shown in 1207 to 1210, till one or more physical links reach full (i.e., 100%) utilization.

Physical links whose capacities are fully allocated are called critical links. Thus the net effect of steps 1207 to 1210 is to recursively identify critical links. While the critical link identification procedure is likely to identify only one physical link as a critical link at each stage of the iterative procedure, the algorithm, as implemented, is equally capable of identifying and processing more than one critical link at a given time.

In one embodiment of the present invention, the identification of critical links is done by incrementing, at 1207, the current Entropy Rate Function estimate by a fixed amount using a functional expression for the Entropy Rate Function that depends on the offered traffic model. Examples of such expressions can be found in Equation (EQN. 15) for homogenous Poissonian traffic, Equation (EQN. 21) for multi-class Poissonian traffic and Equations (EQN. 23 & EQN. 25) for normally-distributed traffic. It should be noted that an increment to the Entropy Rate Function Estimate may, occasionally be negative. This can happen if the dimensioning algorithm overshoots the optimal value and over allocates capacities.

The values of the shift parameter, s, are calculated at 1208 for each of the VPs in the dimensioning set. It should be noted that the shift parameter value represents the slope of the entropy-capacity graph of FIG. 10, for the corresponding VP. The incremental capacities to be allocated to the VPs in the dimensioning set is calculated at 1209 using the incremented values of the Entropy Rate Function. Steps 1207 to 1209 may be performed in a different sequence than that shown in the FIG. 12 based upon implementation considerations.

The capacities allocated to various VPs are then cumulated at 1210 for each of the physical links and compared in 1211 against the total capacity of that physical link. If the unallocated physical capacity of a link falls below a preset limit, then the link is determined to be a critical link.

If the comparison results in the identification of a physical link as a critical link, then the computation proceeds to 1212. If no physical link is found to be a critical link, then steps 1207 to 1210 or their equivalents are iteratively repeated until a critical link is found. In rare situations involving offered traffic models having a truncated right tail, it is possible that the iterative procedure is sometimes unable to identify any more critical links. In such situations, the calculation automatically terminates when the Entropy Rate Function reaches its maximum value, I_(MAX), as specified at 1204.

After identifying a critical link, the VP dimensioning algorithm outputs the results and reformulates the problem as shown at 1212 to 1213. Every time one or more physical links are identified at 1211 to be critical links, the VP dimensioning algorithm proceeds to generate output, at 1212, detailing the currently allocated capacities on each of the VPs traversing a critical link. VPs traversing a critical link are removed from the dimensioning set at 1213. If no VPs remain to be dimensioned, then the dimensioning algorithm terminates at 1216.

If one or more VPs remain to be dimensioned, then the dimensioning set is redefined at 1215 to include only such VPs. Since the VPs traversing critical links were removed from the dimensioning set, and since these removed VPs used up part of the physical link capacities, the dimensioning task reduces to the distribution of the unallocated physical link capacities over the VPs still remaining in the dimensioning set. This is done at 1215 by reducing the available capacities of the various physical links by the amounts allocated to removed VPs that correspond to the critical links last found at 1211. In an alternative embodiment of the invention, the same effect is achieved by freezing the values of the Entropy Rate Function for VPs that have been eliminated from the VP dimensioning set. Since the capacities allocated to removed VPs were generated at 1212, this computation can be easily performed. After the reformulation of the problem at 1215, the algorithm loops back to 1207, and as before, the Entropy Rate Function is incremented by a fixed amount for all the VPs still remaining in the dimensioning set.

Although a preferred embodiment of the method and apparatus of the present invention has been illustrated in the accompanying drawings and described in the foregoing detailed description, it is to be understood that the invention is not limited to the embodiment disclosed, but is capable of numerous rearrangements, modifications and substitutions without departing from the spirit of the invention as set forth and defined by the following claims. 

What is claimed is:
 1. A method for dimensioning virtual paths defined within a telecommunications network carrying general traffic over a plurality of interconnected links whose transmission capacities are limited, said method comprising the steps of: modeling the load on each virtual path of said telecommunications network using a computer configured to implement an entropy rate function; and using said entropy rate function as a blocking measure to operatively cause said computer to dynamically control load balancing of said general traffic over said interconnected links.
 2. The method of claim 1, wherein said step of using said entropy rate function further includes determining actual traffic measurements of at least a portion of said general traffic over certain ones of said interconnected links.
 3. The method of claim 1 wherein said step of using said entropy rate function further includes determining idealized characteristics for said general traffic based on at least one probabilistic function selected from a group comprising a homogenous Poissonian distribution, a multi-class Poissonian distribution, a Normal distribution, and a Binomial distribution.
 4. A method for dimensioning virtual paths on a constrained telecommunications network having a plurality of physical links interconnecting a plurality of nodes, said method comprising the steps of: using a computer to map at least a portion of said plurality of physical links into at least two virtual paths, each of said virtual paths providing an individually switchable connection between a different pair of nodes in the telecommunications network; specifying offered traffic for each of said virtual paths; identifying a transmission capacity constraint for at least a portion of said physical links within said telecommunications network that are available for inclusion in said virtual paths; using said computer to model offered traffic on at least a selected portion of said physical links available for inclusion in said virtual paths based at least in part on an entropy blocking measure function; and selectively allocating capacities to each of said virtual paths, subject to said transmission capacity constraints, based on said entropy blocking measure function such that the blocking probabilities for at least said at least two virtual paths are substantially uniform.
 5. A method for dimensioning a constrained telecommunications network having a plurality of physical links interconnecting a plurality of nodes, said method comprising: mapping a plurality of physical links into one or more virtual paths, each of said virtual paths providing an individually switchable connection between a pair of nodes of the telecommunications network; grouping a selected plurality of said virtual paths into a dimensioning set; allocating initial transmission capacity to each virtual path at an initial blocking value; selecting at least one subsequent blocking value that is lower than a corresponding initial blocking value by: calculating blocking on each of said virtual paths traversing a single physical link; allocating the available transmission capacity amongst virtual paths traversing a single physical link in response to variations in the blocking amongst different virtual paths until a physical link is identified as having no unallocated capacity; removing all virtual paths traversing said identified physical link from said dimensioning set; reducing the transmission capacity allocated to each of the physical links by an amount previously allocated to said removed virtual paths; and iteratively repeating said calculating, allocating, removing and reducing steps until no virtual paths are left in the dimensioning set.
 6. A computer implemented system for dimensioning virtual paths defined on a telecommunications network carrying general traffic, said network having a plurality of interconnected links whose transmission capacities are limited, said system comprising: means for choosing an appropriate entropy rate function to model the load on each virtual path of said telecommunications network; means for selecting a solution algorithm using the entropy rate function as a blocking measure that is operative to solve a load balancing problem for said general traffic; and means for controlling load distribution via a computer using said load balancing algorithm incorporating said entropy rate function to produce a load distribution on said virtual paths that is substantially uniform.
 7. The virtual path dimensioning system of claim 6, wherein said entropy rate function is determined by traffic measurements.
 8. The system of claim 6 wherein said entropy rate function further determines idealized characteristics for said general traffic based on at least one probabilistic function selected from a group comprising a homogenous Poissonian distribution, a multi-class Poissonian distribution, a Normal distribution, and a Binomial distribution.
 9. A system for dimensioning virtual paths on a constrained telecommunications network having a plurality of physical links interconnecting a plurality of nodes, said system comprising: means for mapping said plurality of physical links into one or more virtual paths, each of said virtual paths providing an individually switchable connection between a pair of nodes in the telecommunications network; means for specifying offered traffic for each of said virtual paths; means for specifying a transmission capacity constraint for each physical link of the telecommunications network, including each physical link in each of said virtual paths; means for modeling offered traffic over each of said physical links, using an entropy blocking measure; and at least one computer configured to allocate capacities to each of said plurality of virtual paths based on said modeling and subject to said link transmission capacity constraints, such that the blocking probabilities on the at least a portion of said different virtual paths are significantly uniform.
 10. A computer implemented method for establishing virtual paths within a network configured to carry general traffic, said network having a plurality of interconnected links whose transmission capacities are inherently limited, said method comprising: choosing an appropriate entropy rate function to model the load on each virtual path of said telecommunications network, said entropy rate function being determined by idealizing the characteristics of offered traffic on a telecommunications network, and wherein the idealized entropy blocking measure used to model the offered traffic is said Entropy Rate Function, I_(X)(C), said Entropy Rate Function being calculable as the approximation of the negative logarithm of the probability that an arbitrarily distributed random variable, X, is greater than or equal to a preselected value, C, and said Entropy Rate Function additionally being a convex function obtaining its minimum value of zero at a mean of the distribution; selecting a solution algorithm using said entropy rate function as a blocking measure that is operative to solve a load balancing problem for said general traffic; and performing computations on a computing system using said load balancing algorithm incorporating said entropy rate function to produce a load distribution on said virtual paths that is substantially uniform for at least a selected portion of said virtual paths.
 11. A computer implemented method of dimensioning virtual paths defined on a telecommunications network carrying general traffic, said network having a plurality of interconnected links whose transmission capacities are inherently limited, said dimensioning method comprising the steps of: choosing an appropriate entropy rate function to model the load on each virtual path of said telecommunications network; selecting a solution algorithm using said entropy rate function as a blocking measure that is operative to solve a load balancing problem for said general traffic, wherein said solution algorithm uses said entropy rate function as a blocking measure and is further configured to act as a push down algorithm that assembles the virtual paths to be dimensioned into a dimensioning set, calculates the blocking on each virtual path at every network link traversed by said virtual path using said entropy rate function, identifies the virtual path having the highest blocking on each network link, and allocates additional capacity to the identified virtual path without violating the network resource constraints until it no longer has the highest blocking; and performing computations on a computing system using said load balancing algorithm incorporating said entropy rate function to produce a load distribution on said virtual paths that is substantially uniform.
 12. A computer implemented method of dimensioning virtual paths defined on a telecommunications network carrying general traffic, said network having a plurality of interconnected links whose transmission capacities are inherently limited, said method comprising: choosing an entropy rate function to model a load associated with each of a plurality of virtual paths within said telecommunications network; using said entropy rate function as a blocking measure that is operative to solve a load balancing problem for said general traffic, wherein said algorithm using the entropy rate function as a blocking measure is further configured to perform the following iterative steps: grouping all virtual paths into a dimensioning set; specifying the transmission capacity of each physical link; choosing a relatively large initial blocking value for each virtual path; selecting an error bound to evaluate the convergence of a critical link identifying algorithm; determining the blocking value for each virtual path using an entropy blocking measure; iteratively identifying the virtual path having the largest blocking value as long as the physical link has available capacity; reducing the blocking of the virtual path having the largest blocking value by reallocating transmission capacities amongst said virtual paths as long as no physical link reaches full utilization and the blocking reductions of the maximally-obstructed virtual path is greater than a preselected error bound; identifying physical links lacking allocable capacity as critical links; eliminating all virtual paths spanning critical links from the dimensioning set; and selectively re-adjusting the transmission capacities of the remaining physical links to reflect the capacities allocated to virtual paths most recently eliminated from the dimensioning set; and performing computations on a computing system using said load balancing algorithm incorporating said entropy rate function to produce a load distribution on said virtual paths that is substantially uniform.
 13. A computer implemented system for dimensioning virtual paths defined on a telecommunications network carrying general traffic, said network having a plurality of interconnected links whose transmission capacities are limited, said system comprising: means for choosing an appropriate entropy rate function to model the load on each virtual path of said telecommunications network; means for selecting a solution algorithm using the entropy rate function, and wherein said entropy rate function is determined by idealizing the characteristics of offered traffic on a telecommunications network wherein the idealized entropy blocking measure used to model the offered traffic is the Entropy Rate Function, I_(X)(C), said Entropy Rate Function being calculable as the approximation of the negative logarithm of the probability that an arbitrarily distributed random variable, X, is greater than or equal to a preselected value, C, and said Entropy Rate Function additionally being a convex function obtaining its minimum value of zero at the mean of the distribution; and means for performing computations on a computing system using said load balancing algorithm incorporating said entropy rate function to produce a load distribution on said virtual paths that is substantially uniform.
 14. A computer system for dimensioning virtual paths defined on a telecommunications network carrying general traffic, said network having a plurality of interconnected links whose transmission capacities are limited, said computer system comprising: at least one processor configured to operatively respond to computer implemented instructions associated with an entropy rate function, a solution algorithm, and a load balancing algorithm, wherein: said entropy rate function is configured to model the load on each virtual path of said telecommunications network; said solution algorithm is configured to use the entropy rate function as a blocking measure in solving a load balancing problem for said general traffic by assembling the virtual paths to be dimensioned into a dimensioning set, calculating a blocking on each virtual path at every network link traversed by said virtual path using said entropy rate function, identifying the virtual path having a largest blocking on each network link, and selectively allocating additional capacity to the identified virtual path having the largest blocking to reduce the blocking associated with the identified virtual path; and said load balancing algorithm is configured to incorporate said entropy rate function to produce a load distribution on each of said virtual paths that is substantially uniform.
 15. A system for use in dimensioning virtual paths defined within a telecommunications network for carrying general traffic, said network having a plurality of interconnected links whose transmission capacities are limited, said system comprising: means for choosing an appropriate entropy rate function to model the load on each virtual path of said telecommunications network; means for selecting a solution algorithm that uses said entropy rate function as a blocking measure in solving load balancing problems associated with said general traffic, said solution algorithm iteratively implementing a push down function using: means for grouping all virtual paths into a dimensioning set; means for specifying the transmission capacity of each physical link; means for choosing a relatively large initial blocking value for each virtual path; means for selecting an error bound to evaluate the convergence of a critical link identifying algorithm; means for determining the blocking value for each virtual path using an entropy blocking measure; means for iteratively identifying the virtual path having the largest blocking value as long as the physical link has available capacity; means for reducing the blocking of the virtual path having the largest blocking value by reallocating transmission capacities amongst said virtual paths as long as no physical link reaches full utilization and the blocking reductions of the maximally-obstructed virtual path is greater than a preselected error bound; means for identifying physical links lacking allocable capacity as critical links; means for eliminating all virtual paths spanning critical links from the dimensioning set; and means for iteratively readjusting the transmission capacities of the remaining physical links to reflect the capacities allocated to virtual paths most recently eliminated from the dimensioning set; and means for controlling a load distribution on each of said virtual paths based on said entropy rate function and said solution algorithm such that the resulting load distributions are substantially uniform. 